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Abstract 

In this paper we analyze several anisotropic bootstrap percolation 
models in three dimensions. We present the order of magnitude for the 
metastability thresholds for a fairly general class of models. In our proofs, 
we use an adaptation of the technique of dimensional reduction. We find 
that the order of the metastability threshold is generally determined by 
the 'easiest growth direction' in the model. In contrast to anisotropic 
bootstrap percolation in two dimensions, in three dimensions the order of 
the metastability threshold for anisotropic bootstrap percolation can be 
equal to that of isotropic bootstrap percolation. 
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1 Introduction 

In bootstrap percolation models on a finite cube [0,L] d £ Z d , in the starting 
configuration every site is occupied with probability p, and empty otherwise, 
independent of all other sites. 

The configuration then evolves according to the bootstrap rule: each site that 
has at least k occupied sites in its neighborhood becomes occupied, and occupied 
sites remain occupied. The rule will repeatedly be applied until no new site will 
become occupied. One is usually interested in the probability that the cube is 
internally spanned, that is, in the final configuration all sites are occupied. We 
will in this paper choose k to be half the number of sites in the neighborhood. 
In ordinary bootstrap percolation the neighborhood of a site consists of all its 
nearest neighbors. In anisotropic bootstrap percolation, however, the size or 
shape of the neighborhood is not equal in every direction. For earlier work on 
bootstrap percolation models, see e.g [2j |5J PH P^l EU l |2T | |27 | 125 1 [25]. 

Bootstrap percolation models and arguments have been applied in a variety 
of settings, from fluid dynamics, magnetic models, the theory of glasses, neural 



Figure 1: Illustration of the neighborhood of a site in the (2, 3, 4) model. The z 
axis is towards the reader. 

networks, the theory of sandpiles to rigidity theory and economics, see e.g. 
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Anisotropic bootstrap percolation models until now have been studied in 
two dimensions. For instance, Gravner and Griffeath [T!5] introduced the model 
where the neighborhood consists of six sites, namely, in the x direction only the 
nearest neighbors, but in the y direction both the nearest and the next-nearest 
neighbors. We will call this the (1,2) model. In this model (where k = 3), an 
occupied rectangle can grow by the bootstrap rule most easily in the y direction. 
A single occupied site at distance 1 or 2 from the square suffices to fill the next 
line segment in this direction, whereas in the x direction, two occupied sites 
with no more than three empty sites in between are needed to fill the next line 
segment. The behaviour of this model is similar to that of the semioriented 
model studied in [TJ1 [37J ; for a bootstrap rule whose anisotropy appears to 
be of qualitatively different type, as its asymptotics follows standard isotropic 
behaviour, see [7]. 

We will study the following kind of anisotropic model (see Figure [lj : the 
neighborhood consists of the a-nearest neighbors in the x direction, the ^-nearest 
neighbors in the y direction and the c-nearest neighbors in the z direction, with 
a < b < c, and we choose k = a + b + c. We will call this the (a, b, c) model. 
The notation (a, b) for two-dimensional models, like the (1,2) model of [IS], is 
similar. 

An important tool that we will use is dimensional reduction |28j . Suppose a 
certain large threedimcnsional rectangular block is occupied. Then all sites in 
a slab (we call the set of sites outside the block that are adjacent to one of the 
faces, a slab), already have some number k' — a, b or c of occupied sites in the 
intersection of their neighborhood with the rectangular block. Thus, once the 
rectangular block is occupied, a slab in direction x will become occupied in the 
(a, 6, c) model, if the slab is (internally) spanned in the (6, c) model. For the 
(a, 6, c) model, we will call the (6, c) model the reduced model in the x direction, 
and likewise define the reduced models in the y and z directions. 

Our main example is the (1,1,2) model. The (1,1,2) model has two hard 
growth directions and one easy one. Namely, in the z direction the reduced 
model is (1, 1), but for the other directions it is (1, 2). Our main result however 
is valid for all (a, b, c) models. 
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2 Terminology 



In this section we present some common terms in bootstrap percolation, which 
we will use throughout the paper in the sense explained below. 

Threshold lengths, sharp thresholds Call P([0, L] d internally spanned) 
the probability that for a bootstrap percolation model on a finite cube [0, L] d , in 
the final configuration all sites are occupied. The typical behavior of bootstrap 
percolation models is that there is a sharp threshold length, this means that 
there is a function /(^) such that, as p — > 0, 

P([0,L] d i.s.) 

F([0,L] d i.s.) 
P([0,L} d i.s.) -> 1 

F([0,L] d i.s.) -> 1 

We will call e /( p } (in d = 2), or e 6 ^ 5 ' (in d = 3) the threshold length L*' l (p). 

Inversely, there is a percolation theshold pth{L) = Cg(L), such that asymp- 
totically for increasing L for all positive e there will be percolation if p > 
(1 + s)g(L) and no percolation if p < (1 — e)g(L), with high probability. The 
first proof of a sharp threshold length was for the isotropic model in dimension 
2 by Holroyd [55]. Since then, several more sharp threshold results have been 
obtained [3 E E3 GJ [23 [23 ■ 



for L < e^-^Cp) (in d = 2), or 

for L < e (in d = 3), and 

for L > e (1+e)/( p ) (in d = 2), or 

for L>e e P (hid = 3). 



Order of threshold length There are many bootstrap percolation models 
for which there is no proof (yet) of a sharp threshold length. For some of those 
less precise results known, namely, lower and upper bounds for L th (p) similar to 
those of ( |2.1[ ), but with different multiplicative constants 7 and T > 7, namely: 

L] d i.s.) ^0 for L < L~(p) = e 7/( p } (in d = 2), or 

L] d i.s.)^0 for L < L~(p) = e^*^ (in d = 3), and 
L] d i.s.) 1 for L > L+(p) = e r/( p } (in d = 2), or 

L] d i.s.) -> 1 for L > L+(p) = ef^^ (in d = 3). 

(2.2) 

In such a case we say that the order of the threshold length is known. 
In this paper, we focus on /(|), that is, orders of the threshold length of the 
following kind: 



For (a, b) models, we are interested in the order of In L (p), 
For (a, b, c) models, we are interested in the order of In In L th (p) 



Indeed, our main result is the determination of the order of In In L th (p) for the 
general (a, 6, c) model. However, since strictly speaking we do not prove that a 
sharp threshold length exists, we will not use this term in our theorem. 
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In the remainder of this paper, we will use subscripts a, b and a, b, c to refer 
to the (a, b) model resp. the (a, b, c) model. Sometimes we will use the term 
^abip) m cases where no sharp threshold result is known, that is, f a ^{l/p) is 
known only up to a constant. We will only use this abusive notation in cases 
where the constant is unimportant, to avoid cumbersome elaborations in terms 
of the lower and upper bounds. 

We conjecture that in fact for all the models we consider, that is, every (a, b) 
model in two dimensions, and every (a, b, c) model in three dimensions, there 
exists a sharp threshold length. 

Critical droplets Often, the proof for an upper bound L + (p) for L th {p) 
involves the notion of a "critical droplet" . This is an occupied connected set of 
sites of a size and shape such that, if all other sites in a finite or infinite lattice 
are independently occupied with probability p, then the droplet will continue 
to grow with high probability. Since the occurrence of such a critical droplet is 
(for fixed p) a local event, if L is large enough, that is to say L increases fast 
enough as p — > 0, then there will with high probability be a critical droplet in 
the volume of linear size L. 

As a simple example, consider isotropic bootstrap percolation in one dimen- 
sion, that is, on the interval [0, L]. Suppose that L increases faster than 1/p as 
p — » 0. In this case, the critical droplet is one occupied site. The probability 
that somewhere on the line there is an occupied site, tends to 1. Also, this 
critical droplet grows with probability 1: once there is at least one occupied 
site, then the whole line is spanned. 

Note that there are multiple possibilities for choosing a critical droplet. The 
smaller the size of the critical droplet one chooses, the tighter upper bound 
for L th (p) one obtains. In this paper, our choice is always a certain occupied 
rectangle (rectangular block) . (We expect that this choice of rectangles is not 
optimal; in studies such as |22] for example, the critical droplet is more subtle: it 
is a certain set of occupied sites such that a rectangle of size O(lfp) is internally 
spanned rather than fully occupied.) 

Suppose the number of occupied sites in the critical droplet is V. Then the 
probability P that a fixed site is in a critical droplet, is of order p v . Then, 
disregarding some corrections due to the finite size of the droplet (which are 
negligible if L » V, and p small): 

P([0,L} d i.s.) > 1- (1- P) Ld « 1 - er PL \ 

where we used that (1 — P)v tends to e" 1 as p —> 0. We can alternatively 
interpret this formula as follows: suppose we have a Poisson point process with 
parameter P. Then the above probability is the probability that there is at 
least one Poisson point in the cube [0, L] d . For this reason, one also calls P 
the density of critical droplets. We see from this expression that P([0, L] d i.s.) 
tends to 1 if L d tends to infinity faster than 1/P. We conclude that for every 
P denoting the density of a certain choice of critical droplets, 

L+(p) = 0(P- 1/Cl ). (2.3) 

Inversely, we can say that if L < L~(p), then with high probability there is no 
critical droplet of any kind, so that for every choice of critical droplet, we have 

P = 0{L-{p))- d . (2.4) 
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(a, 6) order of In L^ b (p) 

(1,1) ^ [2J 

(1,6) with 6>1 Mn 2 i[l2] 

(a, a) i first mentioned in [19] (Section 7), see also [TT| H"3] 

Table 1: All known sharp thresholds for (a,b) models. 



Supercritical size We say that a volume is of supercritical size if the prob- 
ability that it is internally spanned, tends to 1 as p — > 0. For example, the 
volume [0,L] d , with L > L + (p), is of supercritical size. However, a volume of 
supercritical size does not need to be cubic. We will also consider rectangular 
volumes of supercritical size. For example, a rectangle is of supercritical size if 
it consists of at least (L + (p)) d sites, and its shape is such that it is much larger 
than the critical droplet, in every direction. 

3 Main result 

Our main result is that the order of lnlnL" l (p) for the general (a,6, c) model, 
depends only on a and b. We state this result as follows: 

Theorem 3.1. For the (a, 6, c) model, there exist constants "f a ,b.c a- n d ^a.b.c 
such that, as p — > 0, 

LY i-s.) for L < e e , 
L\ i.s.) — > 1 for L > e 



where 

• if a = b then f a ,a{ l /p) =P~ a + o(p~ a ), 



of a < b then / a ,b(l/p) = P a In 2 p + o(p a In 2 p) 



Since the proof of this theorem involves dimensional reduction, we use in- 
formation on / a ,h(l/p), the order of the threshold length of the (a, b) model. In 
several cases a sharp threshold length is known, namely, for the (1, 1) model 
[22] and for the (1,6) model with 6 > 1 [12] . In these cases we can specify 
the constant r ai f, !C . For example, it will turn out that if a = 6 = 1, then the 

2 

constant r x l c can be chosen to be twice the Holroyd constant 2Ch — \- 

Furthermore, from the analysis of Duminil-Copin and Holroyd [111 113]. it 
follows that for the (a, a)-model a sharp threshold result holds, namely that 



there is a constant C such that (pD! holds with f a ,b( k ) = e Cp a+o( - p "\ We 



summarize the current knowledge in Table [T] 

For the general case of the (a, 6) model with 1 < a < 6 it is not known 
whether a sharp threshold exists. For the (a, 6) model with a < 6, we derive the 



following result (see Section 4.2), which is sufficient for our purposes: 
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Claim 3.2. Let a <b, then for the (a, b) model, there exist constants ~f a ,b and 
T a b, such that 



{{0,L} 2 is.) ^ for L< L- b (p) = e^"' 1 " 1 ', 
{[Q.L] 2 Is.) ->• 1 /or L > L+ 6 (p) - e^" ^f. 



(3.1) 



Remark. For the (1,1,1) model, a sharp threshold result is already known 
[6j. One might expect that this result should give us a lower bound for 7 a ,&, C j 
since it seems natural to expect (and indeed we do) that an (a, 6, c) model with 
a larger neighborhood should have more difficulty in growing, and in filling up 
a volume. However, a direct inequality is not that obvious. Indeed, as a pair of 
occupied sites on a line parallel to the y-axis at distance 4 in the (l,2)-model 
can cooperate, whereas they cannot in the (1, l)-model (as then they will never 
belong to the neighborhood of the same site) a configuration-wise ordering is 
excluded. 

Thus, we present our intuition as a series of conjectures: 
Conjecture 3.3. 

• If a < b, then there is a constant C a b such that 

\nL%{p) = C a , b p- a ln 2 (l/p) + o(l/p- a ln 2 (l/p)), 

• Ifb' > b, then C a # > C a , 6 . 

Conjecture 3.4. 

• There is a constant C a ,b,c such that 

\nhiLf Ac {p) = C aAc f a , b (l/p) + o(f atb (l/p)), 



where f a ,b(^/p) is as in Theorem 3.1 



If a 1 >a,b'> b, and d > c then In In b , c , (p) > lnlni" 6 c (p) 



We expect that the first half of Conjecture |3.3| might be derived by extend- 
ing the techniques from 12 . Indeed, by the conjecture holds for a = 1. 
Nevertheness, the full proof is outside the scope of the present paper. 

4 The lower bound 

We first explain the lower bound in the case (a, 6, c) = (1, 1, c) in some detail, 
as an instructive example. 
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Figure 2: Illustration of how the critical droplet grows in the (1, 1, c) model (not 
to scale). The z direction is towards the reader. First the droplet grows only in 
the z direction, then in all directions. 

4.1 The lower bound for the (1,1, c) model 

In this section, we show that for all L(p) larger than the value stated in Theorem 
|3.1[ the probability that the configuration is internally spanned, tends to 1. We 
choose an occupied rectangular block as critical droplet. This rectangular block 
will keep growing if the areas of its faces are large enough. For instance, if the 
area of the face in the z direction is larger than e 2C, «( 1 + e )/p ; then the slab next 
to it is spanned with high probability, because its size is supercritical for the 
(1,1) model. 

Naively, one could suppose that we need a critical droplet with one face 
of supercritical size for the (1, 1) model, and two of supercritical size for the 
(1, c) model. However, we can do better than that, following essentially the 
same line of thought as in jTTj which is similar to that in [2] and which was 
based on an unpublished observation of Roberto Schonmann. There a suitable 
critical droplet for the (1, 2) model was found to be a strip of length ^- In -, and 
width 2, with C a large enough constant. Clearly the width of this strip is not 
supercritical, but the strip grows with large probability into a rectangle of size 
p In p times ^ , which is of supercritical size in both directions. The probability 

• •• — In - — — In 2 — 

that such a critical droplet is at a fixed position, is p p p = e v W . 

For the (1,1, c) model, we choose as critical droplet a rectangular block of 
size N by N by 2, with N such that the droplet will grow sufficiently far in the 
z direction. More specifically, let us define 

2C H 

Pi,i=e- — , (4.1) 

which is the probability that a fixed site is in a critical droplet in the (1,1) 

model. Then a square of size N > e » is supercritical for the (1,1) 
model, because with high probability there will be a critical droplet in it. 

We will choose iV 2 = (Pi,i)~( 1+e ' slightly larger than minimally supercrit- 
ical, because we need that many slabs of that size will get occupied with high 
probability, rather than just one. For our critical droplet to grow into a rectan- 
gular block with all faces of supercritical size, we need that M adjacent slabs 
get occupied with large probability, with M such that MN > (L^ c (p)) 2 . As 
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was proved in [171 (and explained above), there is a constant such that 

Lf 2 (p) < e^~ ln p . In [T5], this result was extended: a sharp threshold length 

of the form L^Jjp) = e » ln * was found for all (l,c) models, where C(c) is 
a constant depending only on c. The probability that a slab contains a critical 
droplet for the (1,1) model is at least 1 — (1 — Pi.i) N , and therefore the prob- 
ability that M slabs will get occupied, is at least (1 — (1 — Pi,i) N ) M (this all 
up to some irrelevant constants, due to the finite size of the droplets, possible 
overlapping of the droplets, etc). 

The logarithm of this probability is approximately 

2 - 2 ° H ' 

-Me~ Pl lN = -Me-< p "> ' = -Mer e p , which is close to for all M = 

o(e e p ). We see that it is possible to choose M large enough so that MN > 
(L\ h c (p)) 2 is satisfied. 

Now we need to find a bound for L such that a cube [0, L] 3 contains a critical 
droplet with high probability. This is the case if L 3 > p- 2N (see flsTs])), We 
work out 

2C H (l+e) 2C H (l+e) t 

, 2^2 |i l p j, +o(-) 

L > p 3 = e 3 p = e , 

which gives us the desired lower bound. 

We remark that it is possible to refine this argument and obtain a smaller 
estimate for N. For instance, if we suppose that M — N k for some k, then we 
find that an N by N by 2 occupied rectangular block with iV 2 ~ -p^ ln -pjj— is 
sufficiently large to act as a critical droplet. However, this refinement leads to 
the same conclusion. 

4.2 The lower bound in the general case 

We start with deriving Claim [X2| that is, we derive some basic estimates for 
^ab(p) an d L^bip)- Since the proof mostly consists of repeating work that has 
already been published, here we just sketch the steps that are needed. Namely, 
considering that 

• an occupied rectangle of size at least Cp~ a x b, with C large enough, 
serves as a critical droplet, as such an occupied rectangle grows with high 
probability into an occupied rectangle of size Cp~ a ln | x p~ b , in the same 
manner as the 2 by ln i critical droplet of [17] (the generalization is 
straightforward) , 

• the probability of a fixed 2 x Cp~ a ln | rectangle to be occupied, is p Cp ln ? = 

-Cp~ a In 2 i 

e F p , 

• therefore, b (p) < e Cp ln 

and we obtained one of the inequalities. 

For the bound in the other direction, again a generalization of the method of 



[17] for the (1,2) model is needed. First, we generalize (2) of [T7] to (5.2) (this 



paper, Section 5.1). Then we repeat the calculation that leads to (5) of [17] . 



replacing p 2 by p" and p by p a , and choosing x — 6*2^ In ^ and y = p b+1 / 2 
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instead of choosing k — ^572 and I — C2 ^ In ^ . This will lead to the other 
inequality. 

Now we get back to the (a,b,c) model in three dimensions. We will show 
that for all L larger than the value stated in Theorem |3.1[ the probability that 
the configuration is internally spanned, tends to 1. Suppose that there is an 
occupied rectangular block in the cube [0,L] 3 , then for each direction, the size 
of its face will determine whether it will grow easily in that direction (also, the 
width of the rectangular block should be large enough; width c is sufficient). 
For an occupied rectangular block to serve as a critical droplet, at least one of 
its faces should be of supercritical size for the reduced model in that direction. 
On the other hand, if all of the faces of an occupied rectangle are of supercritical 
size, then it is certainly a critical droplet. 

We suppose for ease of notation that the easiest growth direction is the z 
direction. Only if a < b and Conjecture |3.3| does not hold, then it is possible that 
the (a, c) model has a lower threshold length than the (a, 6) model. But in that 
case, these two models do have the same order of threshold length, therefore, in 
that case the y and z direction would behave essentially in the same manner. 

We denote 

L+(p) = max{L+ (p),L+ e (p),Lt c (p)}. 



Let P a f, denote the density of critical droplets for the (a, b) model (see (2.3 1). 

We choose as a critical droplet for the (a, b, c) model an occupied N by N by 
2 rectangular block, where N > (L"}~ b {p)) 1+e ■ This rectangular block will grow 
easily into a N by N by M rectangular block that is of supercritical size in all 
directions, that is, MN > (L+(p)) 2 . 

The probability that M slabs will get occupied, is at least (1 — (1 — P a ,b) N ) M 
(up to some irrelevant constants). We take the logarithm, which is approxi- 
mately -Me- p «. tJvJ = -Me- p «*. This is close to if M = o{e {L ^ {p)Y ) . By 



the bounds in Claim 3.2 we see that it is possible to satisfy both M > L + (p) 
and A/ = o(e (i ^ (p)r ). 

The probability that a fixed N by N by 2 rectangular block is occupied, is 
p 2N . Therefore, p 2N describes the density of critical droplets in the (a, b, c) 



model, and we obtain (see (2.3)) 

L 3 > (1/ P ) 2N2 . 

We work out 

L > (l/p)§( i+ W) 2+2e = e e«"-..(Vrt(i + .) + .(/..»(t/»)) > (42) 
where C is a suitable constant. This gives the lower bound for L in Theorem 



5 The upper bound 

In this section, we show that for all L smaller than the value L~ bc (p) stated 
in Theorem |3.1[ the probability that the configuration is internally spanned, 
tends to 0. In the case that L < (1/p) 1 / 4 , we simply estimate P([0,L] 3 i.s.) < 
p (p- 1/4)3 _ p i/4 j w hi c h tends to if p 0. In the case that L > (l/p) 1/4 
however, we need to do a lot more work. 
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Figure 3: An example to explain weakly spanned sites. Consider the (2, 3) 
model, in a 6 by 6 rectangle. On the left is the initial configuration (occupied 
sites are marked o), in the middle the final configuration (spanned sites are 
marked s), and on the right we have added weakly spanned sites (marked w). 
We see that the rectangle is weakly crossed. 

We will apply the method introduced in [5] for the case of threedimensional 
isotropic bootstrap percolation, later generalized to arbitrary dimension in [S], 
and again presented in |23) . where some minor inaccuracies in the original pre- 
sentation are corrected. We modify the method to be applicable in the case of 
anisotropic bootstrap percolation. 

A necessary condition for a cube [0, L] 3 to be internally spanned in the 
(a, 6, c) model is that it is internally crossed, that is, in the final set of occupied 
sites, there is a path of occupied sites connecting two opposite faces of the cube, 
in every direction. To find a suitable bound for the probability of [0, L] 3 to 
be internally crossed, we will enhance the configuration, that is, occupy more 
sites. Throughout the remainder, we will suppose for ease of notation that the 
z direction is the easiest growth direction. Again, even though if Conjecture 
|3.3| does not hold then the easiest growth direction might be the y direction, 
this possible exception does not change our proof, since in that case the y and 
z direction behave essentially in the same manner. 

5.1 Some definitions and lemmas 

Definition 5.1. In weakly enhancing the initial configuration in a connected 
set of sites p, where R = [0, r x ] x [0, r y ] x [0, r z ] is the smallest rectangular block 
covering p, we do the following: 

1. we evolve according to the (a, b, c) bootstrap rule, only taking into account 
occupied sites in p. 

2. We occupy sites that are weakly spanned. We find weakly spanned sites as 
follows: for every site that is occupied, but was initially empty, we draw 
lines from this site to every site in its neighborhood that is occupied in the 
final configuration. A weakly spanned site is an empty site that is crossed 
by (at least one) such a line. We give an example in Figure [3j 

We say a rectangular block is weakly crossed if it is internally crossed in 
every direction after weakly enhancing. 

The following lemma is somewhat similar to Lemma 1 from [3]. Related 
arguments continue to reappear in various forms in many papers on bootstrap 
percolation. 
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Lemma 5.2. For the (a,b,c) model, suppose that the cube [0, I] 3 is internally 
spanned. Then there are constants k and A that only depend on a, b and c, 
such that for every k G [1, there is a rectangular block with longest side in 

[k, nk + A] that is weakly crossed. 

Proof. First, we need some further notation. We say a region is a set of occupied 
sites such that the set is connected after weakly enhancing it. Therefore, if p 
is a region then the smallest rectangular block that covers it, is weakly crossed. 
We say that an empty site x has a set of influencing regions if there is a set 
of regions such that if we empty the entire configuration except these regions 
then x is spanned, but if we furthermore empty any one of these regions then 
x is not spanned. There may be several possibilities for the set of influencing 
regions. The cardinality of these sets is always between 1 and a + b + c, and 
the distance between any two regions in one of these sets is at most 2c + 1. We 
say that the influencing rectangular block is the smallest rectangular block such 
that it covers the set of influencing regions. 

Next, we present a special order of occupying sites according to the bootstrap 
rule. After step i, we will have a set of regions Cj. We start with the set C\ of 
regions with diameter 1, that consists of all single occupied sites in the initial 
configuration. Since we suppose that [0,l] 3 is internally spanned, C\ is not 
empty. Then we iterate in the following manner: 

• If there are two regions such that their union is again a region, then d + \ 
consists of all the regions in Ci with these two regions removed, and with 
a new region added that is the union of these two regions. If there are 
several such pairs of regions, then we make an arbitrary choice. 

• If not, then we choose an arbitrary site x that has a set of influencing 
regions. If it has several such sets, then we make an arbitrary choice. Call 
this set Si. 

• We suppose the entire configuration is empty except the regions in Si. 

• We occupy x. Either we now have a single occupied connected component, 
or we would have one after weakly enhancing. Therefore, the union of all 
regions in Si, plus x, forms a region. 

• Ci + i consists of all the regions in d with the regions in Si removed, and 
with a new region added that consists of the union of the regions in Si, 
and x. 

Now the following holds: If the maximal diameter of all regions in Ci is Di, 
then the maximal diameter of all regions in C^+i is in [Di, (a + b+c)(Di + 2c+l)], 
and since we supposed that [0, 1] 3 is internally spanned, the iteration does not 
end until the maximal diameter is I. Therefore, if we choose k = a + b + c, and 
A = (2c + l)(a + b + c), the proof is complete. □ 

Definition 5.3. In enhancing the initial configuration in a cube [0, 1] 3 , we do 
the following: 

1. We divide the cube in l/s disjoint slices [0, l) 2 [i, i + s], with i = . . . l/s—l. 
We suppose for ease of notation that I is a multiple of s. We choose s 



11 



as the minimal value such that if an empty site x is in slice i, and all 
other slices are fully occupied, then site x still needs a + b occupied sites 
in the intersection of its neighborhood with slice i, to become occupied. 
For instance, in J3 \3j the slices have width 2, whereas for the modified 
model in [23] the slices can have width 1. For the (a, &, c) model, we need 
that a + b + c — (2c + 1 — s) = a + b, this implies that we have to choose 
s = c + 1. We say that a set of sites (x, y, i), (x,y,i + 1), . . . , (x,y, i + s) 
is a minicolumn in slice i. For every minicolumn in every slice, if at 
least one site in the minicolumn is occupied, then we occupy all sites in 
the minicolumn. We say that q is the probability for a minicolumn to be 
occupied, and note that p < q < sp, since in the initial configuration, 
minicolumns are independent. 

We evolve the minicolumn configuration of every slice by the twodimen- 
sional bootstrap rule, that is, we view minicolumns as sites in the (a, b) 
model with q as parameter, and iteratively occupy those that have at least 
half of their neighborhood occupied. By our choice of s, we are not occu- 
pying sites that would not get occupied in the threedimensional model. In 
fact, notice that by now we have occupied at least every site that would also 
get occupied with the threedimensional bootstrap rule, and also a subset of 
the weakly spanned sites, namely the ones that are crossed by a line in the 
z-direction. 

We fully occupy every slice where there is an occupied component of at 
least some size S. We call such a slice flooded. We will choose S such 
that it resembles a critical droplet for the (a, b) model with q as parameter. 

We occupy weakly spanned minicolumns. That is, we say a minicolumn 
is weakly spanned if it contains a weakly spanned site. This step is not 



essential for arriving at (5.6), but by including this step we ensure that an 



enhanced configuration dominates a weakly enhanced configuration, which 



is necessary to arrive at (5. 



5. We fully occupy slice and slice L/s — 1. 

We say a rectangular block is e-crossed if it is internally crossed in the z- 
direction after enhancing. We denote Pi(e-crossed) for the probability that the 
cube [0, 1] 3 is e-crossed. 

The following lemma bears resemblance to the derivation of (3.30) of [5], 
and of (23) of [23], but we derive the bound for the (a, b) model with arbitrary 
a and b. In the case (a, b) = (1, 2), one may choose x c < | In ^ and y c < p~ 3 / 2 
as in |17] . but note that this choice for y c is not essential: the lemma is valid 
for any y c < p~ 2+t . 

Lemma 5.4. For the (a, b) model and p small enough, the expected size x °f 
the occupied component of the origin, given that x < x c and y < y c , with 
%c < He < p~ b+e for a fixed e > 0, is at most ^Jp. This bound does not change 
when weakly spanned sites are occupied. 

Proof. We write 

Xc Vc 

X < xyP( diameter = max{x, y}). 

x—l y—1 
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For arbitrary r, we have 

X < r 2 P(0 < diameter < r) + x c y c P(r < diameter < max{a; c , y c })- 

We have that if the diameter of the occupied component of the origin is at 
most r, then there must be at least one occupied site among at most r 2 sites, 
therefore, P(0 < diameter < r) < r 2 p. If the diameter is in [r, max{cc c , y c }], then 
we estimate P(r < diameter < max{x c , y c }) < P( diameter < max{a; c , y c }). To 
estimate this probability, we generalize (2) of [T7] . We first quote this equation, 
which was used for the case (a, b) — (1,2): 

P(a fixed x,y rectangle is internally spanned) < min{(l— (1— p 2 ) y ) x , (1— (1— p) x ) y }, 

(5.1) 

where p > p and p > p are of the order p. In words, the meaning of this 
expression is that a necessary condition for a rectangle to be internally spanned, 
is that there is no row or column in which not even one site is spanned, even if 
the neighboring row or column is fully occupied. In the y direction, one needs 
two occupied sites, close enough to each other, for a third one to be spanned. 
The probability for this to occur at a fixed position, is of order p 2 . In the x 
direction however, only one occupied site suffices. Therefore, the generalization 
to arbitrary a, b is straightforward: 

P(a fixed x,y rectangle is internally spanned) < min{(l— (1— p b ) v ) x , (1— (1— p a ) x ) y }- 

(5.2) 

However, we need a bound that is also valid when weakly spanned sites are 
occupied as well. Note that for a site to be weakly spanned, there has to be a 
spanned neighbor in a nearby row/column, that is, at most a resp. b columns 
resp. rows away. Therefore, we adapt the bound as follows: 

P(a fixed x, y rectangle is internally spanned) < 

min{(l - (1 - p&)(*»+l)v)* (1 _ (l-p a )( 2b+1 ^)y}. 

(5.3) 

We estimate 

min{(l - (1 - p6)(2a+l)»)x ) (1 - (1 - 0a)(»+l)x)v} < (1 - (1 - f)(2a+l)vy 

< (2(2a+l)yp b ) x 

< (2(2a+l)y c p b ) r . 

(5.4) 

Now we insert that y c < p~~ b+e , so that 

P( diameter < ma,x{x c ,y c }) < (2 (2a + l)p e ) r ■ 

Therefore, we find 

X<r i p + x c y c (2(2a + \)f) r < r 4 p + p- 2b+2e (2(2a + \)f) r . 

We choose r such that —26 + (2 + r)e > 1/2. With this choice of r and p small 
enough, x is l ess than ^/p, as needed. □ 
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Lemma 5.5. Let P a ^ = (L b (p)) . Let P/(e- crossed) and s be defined as in 
Definition \5.3\ Then 



'2l 2 {2^Ip) l ' s - 2 «/(l/p) 1/4 < I <p~ 



b+e 



(e-crossed) << / x \i/s 



Al 



{l 3 (sp)^ b P a<b )) 3 ifl>p- b+ \ 



Proof. Since we consider the cube [0, /] 3 to be e-crossed, we denote according 
probabilities with a subscript I. 

We condition on the number of flooded slices: 

l/s-2 

P; (e-crossed) = P ; (e-crossed|m slices flooded)P(m slices flooded). 

m=Q 

Note that the first and last slice are fully occupied, but they are not flooded 
slices. By independence of the slices, we have 

P/(m slices flooded) = (P/(a slice is flooded))" 1 , 

where the probability for a slice to be flooded is the probability that there is an 
occupied component of at least size S for the (a, b) model in the slice. We now 
condition on the positions of the flooded slices: 

P;(e-crossed|ra slices flooded) = P; (e-crossed | slices «i, . . . , i m flooded). 

i\ .. . .,i m 

By £i j, for j > i + we denote the event that there exists a crossing from slice 
i + 1 to slice j — 1. Then 

m 

P; (e-crossed | slices h,...,i m flooded) = Pj(£j^ +1 ). 

j=o 

The slices ij + 1, . . . ,ij+i — 1 are not flooded, therefore the largest occupied 
connected component is strictly smaller than S. If there is a crossing, then 
there is a self-avoiding path from slice ij to slice ij+i that visits a connected 
component in each slice in between, and possibly in some slices more than once 
through several disjoint connected components, like skipping between slices from 
one stepping stone to another. We choose the path such that the total number 
H of connected components visited is minimal. Let Ij = ij+\ — ij — 2. We sum 
over possible values of H : 

H£y,i s+ J < J2 P i(^ J+ i in h ste P s )- 

h>lj 

Note that H can take only values = Ij + 2k, with k = 0, 1, Since for each 

larger value of k the path contains an extra back- and forwardstep, for each k 
the number of possibilities for the order of slices visited is (Ij) . We bound this 
by 2 hk , which is the number we get if we suppose that for every skip there are 
two possibilities, disregarding the fact that the path needs to lead to slice ij+i- 
We need the probability of H connected components: after each stepping stone, 
we need a new stepping stone at the right position for the path to continue. 
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That is, if a stepping stone is contained in a rectangle of x times y, then there 
are at most xy choices for the site from where we skip to the next slice. Since 
slice ij is fully occupied, there are I 2 choices for the site from which we make the 
first skip. We also use that a stepping stone is smaller than x c y c , the dimensions 
of S. We obtain 

h k 

- ^ E 2 ' ifc II E x iVi^l (occupied component has size 

fe>l i=l Xi<x c ,yi<y c 

The factor ^2 X <X y . <y x^y^i (occupied component has size £,2/,) does not de- 
pend on i, and is the expected size x °f the connected occupied component of a 
fixed site, given that it is smaller than x c y c . We will choose S such that we can 



bound this expected size by yfq. By Lemma 5.4 we can choose x c < y c < p~ b+t 
for some fixed e > 0, so that ds = p~ , where ds is the diameter of S. We 
work out 

fc>l i=l fc>l 

where the last bound is valid for q small enough so that 2^/q < 1/2. Therefore, 

F l (£ ijtij+1 )<2f(2^, (5.5) 

If 5 is a critical droplet, then we can estimate P;(a slice is flooded) < {dsl) 2 P a .b 
if I > c?5, and otherwise, since P a ,b is an upper bound for the probability of 
any choice of critical droplet (see (2.4)). However, to be sure that S is a critical 
droplet we need ds > p~ b , but we chose ds < p~ b+e . The probability that S is 
occupied is therefore at most a factor p~ 2bt larger than P a ,b, so we estimate 

P z (a slice is flooded) < (d s lfp~ 2h( P a ,b- 

Putting everything together, we get that if I > p~ b+e , then 

l/s-2 m 

P,(e-crossed) < £ ((d s l) 2 p- 2be P a , b ) m ]T 4Z 2 JJ g^/ 2 . 

m—0 ii,...,i m j— 

We use that the number of possibilities for ii, . . . ,i m is at most / m , and that 
+ I ■ ■ ■ + l m = I - m, so that njLo 1 h/2 = q 1/2{l ~ m) - Then 



(e-crossed) < 4?V /2 £ ( (d s /) V 26e P a ,6 



m=0 



so that 

Pz(e-crossed) < 4Z 2 ^/ 2 ( -^-{d s l) 2 p~ 2b€ P c 
We insert ds < f» 6 and q < sp and simplify, so that we finally get, for Z > p~ b+e , 
P/(e-crossed) < 42 2 (z 3 (sp)^p 26(1 ~ e) P , 6 )Y * . (5.6) 
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In the case that I < p b+e , we have that P/(a slice is flooded) = 0. Therefore 
we use (5.5) with ij = 1 and ij+i = l/s — 1, and obtain 

P/ (e-crossed) < 2l 2 {2^p) l/s - 2 , (5.7) 

which tends to for all (l/p) 1/4 < I < p~ b+£ . 

□ 



5.2 Proof for the upper bound 

Proof of Theorem\3.1\ upper bound for L 



We use Lemma 



5.2 



to estimate the probability of the cube [0,L] 3 to be 



internally spanned as follows: 

),L] 3 i.s.) < L 3 min ((k— 1)A; + A) max P; (e-crossed). (5-8) 



To arrive at this expression, we have reasoned as follows: If a rectangular 
block is weakly crossed, then the smallest cube covering this rectangular block 
has a crossing in some direction, after weakly enhancing. The probability for 
this is bounded by the probability that there is a crossing in the easiest direction 
after weakly enhancing. We will insert our upper bound for Pj (e-crossed) from 
Lemma |5.5| This bound was derived for the enhanced configuration, which 
dominates the weakly enhanced configuration. 

We consider three ranges for L. 

If L < p~ b+e , then we estimate 

P([0,L] 3 i.s.) < P L (e-crossed), 



which tends to by (5.7 1. 

Next, if L > p~ b+e then we combine equations (5.8) and (5.6), and get 

P([0,L] 3 i.s.) < L 3 min ((/t-l)fc+A) max Al 2 (z 3 (sp)^p 2b(1 ~ £) P aib )Y /S 

Kfe< p — * k<l<Kk-\-X V ' / 

(5.9) 

We will use that for p small enough, {sp)^s~ p 2b ( 1 ~ e ) < l. 
In the case that p~ b+<i < L < (L"^)) 1 / 4 , we choose k = ^77^, then (k — 
l)k + A < L. This gives 



([0, L] 3 i.s.) < L 4 max Al 2 (l 3 P a b )) 

i=A<Z<L 



l/s 



We estimate further 

13 ■ r ,\ / a ™„„ r6 /r3n _\\ l / s 



P([0,L] 3 i.s.) < 4 max L 6 (L 3 P ab )) 

1^±<1<L 

< 4 max (L 9 P a , fc )) 



-<1<L 



Recall that (L^,») 1/3 = P a ,b)- H L < (i^ b (p)) 1/4 , then (L 9 P a , b ) tends to 
as p — > 0, so that this entire expression tends to 0. 
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Finally, we consider the case L > (L a ^(p)) 1 ^ 4 - We now choose k = ^"'''^ 

in @, so that (L ^» [p) ^'^ x < l < (.^(p)) 1 /*; this means that I = 0((£~ b (ij)) 1/4 ). 
We also have (k - l)fc + A < (£~ b ) 1/4 , so that: 

L] 3 i.s.) < L 3 max ; (4i- b (p)) 3 / 4 (Z 3 P a , b )' /s 
<4L 3 ma Xi ((^(p))" 3 / 2 )' 7 * . 



We write L~ 6 (p) = e /(1/p) . Since i = 0((£~ b (p)) 1/4 ), there is a constant C 

_ Of a (,(l/p) 

such that the right factor of this expression tends to as e e . Therefore, 

if L 3 tends to oo more slowly, then the entire expression tends to 0. This gives 
the required upper bound for L. 

□ 
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